Advanced tracking systems design and analysis

The results of an assessment of several types of high-accuracy tracking systems proposed to track the spacecraft in the National Aeronautics and Space Administration (NASA) Advanced Tracking and Data Relay Satellite System (ATDRSS) are summarized. Tracking systems based on the use of interferometry and ranging are investigated. For each system, the top-level system design and operations concept are provided. A comparative system assessment is presented in terms of orbit determination performance, ATDRSS impacts, life-cycle cost, and technological risk

Validation and assessment of variant calling pipelines for next-generation sequencing

Background: The processing and analysis of the large scale data generated by next-generation sequencing (NGS) experiments is challenging and is a burgeoning area of new methods development. Several new bioinformatics tools have been developed for calling sequence variants from NGS data. Here, we validate the variant calling of these tools and compare their relative accuracy to determine which data processing pipeline is optimal. Results: We developed a unified pipeline for processing NGS data that encompasses four modules: mapping, filtering, realignment and recalibration, and variant calling. We processed 130 subjects from an ongoing whole exome sequencing study through this pipeline. To evaluate the accuracy of each module, we conducted a series of comparisons between the single nucleotide variant (SNV) calls from the NGS data and either gold-standard Sanger sequencing on a total of 700 variants or array genotyping data on a total of 9,935 single-nucleotide polymorphisms. A head to head comparison showed that Genome Analysis Toolkit (GATK) provided more accurate calls than SAMtools (positive predictive value of 92.55% vs. 80.35%, respectively). Realignment of mapped reads and recalibration of base quality scores before SNV calling proved to be crucial to accurate variant calling. GATK HaplotypeCaller algorithm for variant calling outperformed the UnifiedGenotype algorithm. We also showed a relationship between mapping quality, read depth and allele balance, and SNV call accuracy. However, if best practices are used in data processing, then additional filtering based on these metrics provides little gains and accuracies of >99% are achievable. Conclusions: Our findings will help to determine the best approach for processing NGS data to confidently call variants for downstream analyses. To enable others to implement and replicate our results, all of our codes are freely available at http://metamoodics.org/wes

A Hybrid Likelihood Model for Sequence-Based Disease Association Studies

In the past few years, case-control studies of common diseases have shifted their focus from single genes to whole exomes. New sequencing technologies now routinely detect hundreds of thousands of sequence variants in a single study, many of which are rare or even novel. The limitation of classical single-marker association analysis for rare variants has been a challenge in such studies. A new generation of statistical methods for case-control association studies has been developed to meet this challenge. A common approach to association analysis of rare variants is the burden-style collapsing methods to combine rare variant data within individuals across or within genes. Here, we propose a new hybrid likelihood model that combines a burden test with a test of the position distribution of variants. In extensive simulations and on empirical data from the Dallas Heart Study, the new model demonstrates consistently good power, in particular when applied to a gene set (e.g., multiple candidate genes with shared biological function or pathway), when rare variants cluster in key functional regions of a gene, and when protective variants are present. When applied to data from an ongoing sequencing study of bipolar disorder (191 cases, 107 controls), the model identifies seven gene sets with nominal p-values<0.05, of which one MAPK signaling pathway (KEGG) reaches trend-level significance after correcting for multiple testing. © 2013 Chen et al

CSSs in a sample of B2 radio sources of intermediate strength

We present radio observations of 19 candidate compact steep-spectrum (CSS) objects selected from a well-defined, complete sample of 52 B2 radio sources of intermediate strength. These observations were made with the VLA A-array at 4.835 GHz. The radio structures of the entire sample are summarised and the brightness asymmetries within the compact sources are compared with those of the more extended ones, as well as with those in the 3CRR sample and the CSSs from the B3-VLA sample. About 25 per cent of the CSS sources exhibit large brightness asymmetries, with a flux density ratio for the opposing lobes of $>$5, possibly due to interaction of the jets with infalling material. The corresponding percentage for the larger-sized objects is only about 5 per cent. We also investigate possible dependence of the flux density asymmetry of the lobes on redshift, since this might be affected by more interactions and mergers in the past. No such dependence is found. A few individual objects of interest are discussed in the paper.Comment: 10 pages, 7 figures, 2 tables; accepted for publication in Astronomy and Astrophysic

Are genetic risk factors for psychosis also associated with dimension-specific psychotic experiences in adolescence?

Psychosis has been hypothesised to be a continuously distributed quantitative phenotype and disorders such as schizophrenia and bipolar disorder represent its extreme manifestations. Evidence suggests that common genetic variants play an important role in liability to both schizophrenia and bipolar disorder. Here we tested the hypothesis that these common variants would also influence psychotic experiences measured dimensionally in adolescents in the general population. Our aim was to test whether schizophrenia and bipolar disorder polygenic risk scores (PRS), as well as specific single nucleotide polymorphisms (SNPs) previously identified as risk variants for schizophrenia, were associated with adolescent dimension-specific psychotic experiences. Self-reported Paranoia, Hallucinations, Cognitive Disorganisation, Grandiosity, Anhedonia, and Parent-rated Negative Symptoms, as measured by the Specific Psychotic Experiences Questionnaire (SPEQ), were assessed in a community sample of 2,152 16-year-olds. Polygenic risk scores were calculated using estimates of the log of odds ratios from the Psychiatric Genomics Consortium GWAS stage-1 mega-analysis of schizophrenia and bipolar disorder. The polygenic risk analyses yielded no significant associations between schizophrenia and bipolar disorder PRS and the SPEQ measures. The analyses on the 28 individual SNPs previously associated with schizophrenia found that two SNPs in TCF4 returned a significant association with the SPEQ Paranoia dimension, rs17512836 (p-value=2.57x10-4) and rs9960767 (p-value=6.23x10-4). Replication in an independent sample of 16-year-olds (N=3,427) assessed using the Psychotic-Like Symptoms Questionnaire (PLIKS-Q), a composite measure of multiple positive psychotic experiences, failed to yield significant results. Future research with PRS derived from larger samples, as well as larger adolescent validation samples, would improve the predictive power to test these hypotheses further. The challenges of relating adult clinical diagnostic constructs such as schizophrenia to adolescent psychotic experiences at a genetic level are discussed

Mapping gene associations in human mitochondria using clinical disease phenotypes

Nuclear genes encode most mitochondrial proteins, and their mutations cause diverse and debilitating clinical disorders. To date, 1,200 of these mitochondrial genes have been recorded, while no standardized catalog exists of the associated clinical phenotypes. Such a catalog would be useful to develop methods to analyze human phenotypic data, to determine genotype-phenotype relations among many genes and diseases, and to support the clinical diagnosis of mitochondrial disorders. Here we establish a clinical phenotype catalog of 174 mitochondrial disease genes and study associations of diseases and genes. Phenotypic features such as clinical signs and symptoms were manually annotated from full-text medical articles and classified based on the hierarchical MeSH ontology. This classification of phenotypic features of each gene allowed for the comparison of diseases between different genes. In turn, we were then able to measure the phenotypic associations of disease genes for which we calculated a quantitative value that is based on their shared phenotypic features. The results showed that genes sharing more similar phenotypes have a stronger tendency for functional interactions, proving the usefulness of phenotype similarity values in disease gene network analysis. We then constructed a functional network of mitochondrial genes and discovered a higher connectivity for non-disease than for disease genes, and a tendency of disease genes to interact with each other. Utilizing these differences, we propose 168 candidate genes that resemble the characteristic interaction patterns of mitochondrial disease genes. Through their network associations, the candidates are further prioritized for the study of specific disorders such as optic neuropathies and Parkinson disease. Most mitochondrial disease phenotypes involve several clinical categories including neurologic, metabolic, and gastrointestinal disorders, which might indicate the effects of gene defects within the mitochondrial system. The accompanying knowledgebase (http://www.mitophenome.org/) supports the study of clinical diseases and associated genes

Negative pressure characteristics of an evaporating meniscus at nanoscale

This study aims at understanding the characteristics of negative liquid pressures at the nanoscale using molecular dynamics simulation. A nano-meniscus is formed by placing liquid argon on a platinum wall between two nano-channels filled with the same liquid. Evaporation is simulated in the meniscus by increasing the temperature of the platinum wall for two different cases. Non-evaporating films are obtained at the center of the meniscus. The liquid film in the non-evaporating and adjacent regions is found to be under high absolute negative pressures. Cavitation cannot occur in these regions as the capillary height is smaller than the critical cavitation radius. Factors which determine the critical film thickness for rupture are discussed. Thus, high negative liquid pressures can be stable at the nanoscale, and utilized to create passive pumping devices as well as significantly enhance heat transfer rates

Genome-wide linkage analysis of 972 bipolar pedigrees using single-nucleotide polymorphisms.

Because of the high costs associated with ascertainment of families, most linkage studies of Bipolar I disorder (BPI) have used relatively small samples. Moreover, the genetic information content reported in most studies has been less than 0.6. Although microsatellite markers spaced every 10 cM typically extract most of the genetic information content for larger multiplex families, they can be less informative for smaller pedigrees especially for affected sib pair kindreds. For these reasons we collaborated to pool family resources and carried out higher density genotyping. Approximately 1100 pedigrees of European ancestry were initially selected for study and were genotyped by the Center for Inherited Disease Research using the Illumina Linkage Panel 12 set of 6090 single-nucleotide polymorphisms. Of the ~1100 families, 972 were informative for further analyses, and mean information content was 0.86 after pruning for linkage disequilibrium. The 972 kindreds include 2284 cases of BPI disorder, 498 individuals with bipolar II disorder (BPII) and 702 subjects with recurrent major depression. Three affection status models (ASMs) were considered: ASM1 (BPI and schizoaffective disorder, BP cases (SABP) only), ASM2 (ASM1 cases plus BPII) and ASM3 (ASM2 cases plus recurrent major depression). Both parametric and non-parametric linkage methods were carried out. The strongest findings occurred at 6q21 (non-parametric pairs LOD 3.4 for rs1046943 at 119 cM) and 9q21 (non-parametric pairs logarithm of odds (LOD) 3.4 for rs722642 at 78 cM) using only BPI and schizoaffective (SA), BP cases. Both results met genome-wide significant criteria, although neither was significant after correction for multiple analyses. We also inspected parametric scores for the larger multiplex families to identify possible rare susceptibility loci. In this analysis, we observed 59 parametric LODs of 2 or greater, many of which are likely to be close to maximum possible scores. Although some linkage findings may be false positives, the results could help prioritize the search for rare variants using whole exome or genome sequencing

Transcriptome-based polygenic score links depression-related corticolimbic gene expression changes to sex-specific brain morphology and depression risk

Studies in post-mortem human brain tissue have associated major depressive disorder (MDD) with cortical transcriptomic changes, whose potential in vivo impact remains unexplored. To address this translational gap, we recently developed a transcriptome-based polygenic risk score (T-PRS) based on common functional variants capturing ‘depression-like’ shifts in cortical gene expression. Here, we used a non-clinical sample of young adults (n = 482, Duke Neurogenetics Study: 53% women; aged 19.8 ± 1.2 years) to map T-PRS onto brain morphology measures, including Freesurfer-derived subcortical volume, cortical thickness, surface area, and local gyrification index, as well as broad MDD risk, indexed by self-reported family history of depression. We conducted side-by-side comparisons with a PRS independently derived from a Psychiatric Genomics Consortium (PGC) MDD GWAS (PGC-PRS), and sought to link T-PRS with diagnosis and symptom severity directly in PGC-MDD participants (n = 29,340, 59% women; 12,923 MDD cases, 16,417 controls). T-PRS was associated with smaller amygdala volume in women (t = −3.478, p = 0.001) and lower prefrontal gyrification across sexes. In men, T-PRS was associated with hypergyrification in temporal and occipital regions. Prefrontal hypogyrification mediated a male-specific indirect link between T-PRS and familial depression (b = 0.005, p = 0.029). PGC-PRS was similarly associated with lower amygdala volume and cortical gyrification; however, both effects were male-specific and hypogyrification emerged in distinct parietal and temporo-occipital regions, unassociated with familial depression. In PGC-MDD, T-PRS did not predict diagnosis (OR = 1.007, 95% CI = [0.997–1.018]) but correlated with symptom severity in men (rho = 0.175, p = 7.957 × 10−4) in one cohort (N = 762, 48% men). Depression-like shifts in cortical gene expression have sex-specific effects on brain morphology and may contribute to broad depression vulnerability in men